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Abstract 

Although about 40% of the soft X-ray background emission in 0.4 to 1 keV range has extragalactic 
origins and thus is totally blocked by the Galactic absorption in midplane directions, it decreases at most by 
about 20% in midplane. Suzaku observation of the direction, (£, b) = (235, 0), showed an Ovii emission 
intensity comparable with that of the MBM-12 on cloud Suzaku observation, but revealed a narrow bump 
peaked at ~ 0.9 keV. The latter component is partly filling the decrease of the extragalactic component 
in midplane. The feature can be well represented by a thin thermal emission with a temperature of about 
0.8 keV. Because of the high pressure implied for spatially extended hot gas, the emission is likely a sum of 
unresolved faint sources. We consider a large fraction of the emission originates from faint dM stars. We 
constructed a model spectrum for spatially unresolved dM stars that consistently explains the observed 
spectrum and the surface brightness. The model also suggests that the emission from dM stars decreases 
very rapidly with increasing b, and thus that it cannot compensate entirely the decrease of the extragalactic 
component at 6~2 — 10°. 

Key words: Galaxy: disk — Galaxy: stellar content — X-rays: diffuse background — X-rays: ISM 
— X-rays: stars 
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1. Introduction 

The soft X-ray sky below 1 keV is spatially smooth af- 
ter subtracting the local structures, such as Loop I. In 
high Galactic latitudes, ~ 40 % of the emission is at- 
tributed to emission from faint extragalactic objects, i.e. 
the Cosmic X-ray Background (CXB), in the ROSAT R45 
band ( - 0.44 - 1.0 keV) (McCammon et al. 2002). 
The rest of the emission is considered to consist of emis- 
sion lines from hot gas in the disk and halo of our galaxy 
(McCammon et al. 2002), and from the Heliosphere by 
the solar wind charge exchange (SWCX) process (Cox 
1998; Cravens 2000; Lallement 2004). A small frac- 
tion can arise from intergalactic space. In the Galactic 
midplane, the interstellar X-ray absorption column den- 
sity, Nn, is ~ lO^^cm"^ even in the anti-center direction. 
Therefore, the extragalactic X-ray photons below 1 keV 
are totally blocked. Nevertheless, the R45 band X-ray 
surface brightness decreases only by 20 % or less from 
high Galactic latitude to midplane. This issue has been 
known as the "M band problem" (McCammon & Sanders 
1990; Cox 2005). The M band is the name of a simi- 
lar energy band in the Wisconsin and the Nagoya-Leiden 
rocket programs. Since X-ray photons below 1 keV can 
travel only about 1 kpc in the Galactic disk, there must be 
emission in the midplane within 1 kpc which compensates 
partly the decrease of the extragalactic emission. Nousek 
et al. (1982) and Sanders et al. (1983) suggested hot gas 
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of 3 X 10^ K as the origin, while Rosner et al. (1981) 
pointed out emission from dM stars can contribute ~ 20 
% of the total diffuse emission. Cox (2005) showed that 
if a significant fraction 1/2) of emission originates from 
hot gas in the temperature range of 2.5 x 10^ to 6.3 x 10^ 
K, the hot gas must expand because of its high pressure. 
He suggested young expanding superbubbles or supernova 
remnants evolving in low density region as candidates for 
the emission. 

As an example, the surface brightness averaged over 
rectangular areas of {A£,Ab) = (10°, 2°) along the line of 
£ = 235° is plotted as a function of b in Figure 1. A model 
surface brightness profile consisting of an unabsorbed con- 
stant emission and the CXB absorbed by the average col- 
umn density is plotted together with the observational 
data. For \b\ < 10°, there is 20 x IQ-^c s-^amin'^ of 
excess over the model. This corresponds to about 20% 
of the total diffuse emission at high latitudes. Note that 
the model curve is for the minimum possible absorption: 
if the other 60% of the R45 emission at high latitudes 
is produced by a hot halo, then this too would be largely 
absorbed in the plane, depending on its scale height struc- 
ture. The surface brightness profile suggests asymmetry 
between the profiles for 6 > and 6 < 0, which is more 
pronounced in R4 band. In this paper, however, we con- 
centrate on the excess fiux at b = 0. In section 3, we 
will construct a model which can consistently explain the 
excess at 6 = 0, although we will find the model cannot 
explain the excess in 6 ~ 2 — 10°. 

The origin of the excess midplane emission is not known 
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Fig. 1. ROSAT diffuse X-ray R45 band map (Snowden et al. 1997), and the surface brightness and neutral hydrogen col- 
umn density as functions of b along £ = 235°. The thick white circle in the map indicates the pointing direction of the 
present observation. The surface brightness and the hydrogen column density were averaged over rectangular areas of size 
(A£, Afe) = (10°, 2°). The surface brightnesses in R4, R5, and R45 bands are plotted as step functions. A model surface brightness 
for R45 band which consists of an unabsorbed constant emission (65 X 10~^c s~^amin~^) and the cosmic X-ray background emission 
(10 (E/1 keV)"^ ** photons s~^ cm"'^ str~^ keV~^ ) absorbed by the average column density is shown with a thick curve. There 
exists about 20 X 10~^c s^'^amin"^ of excess over the model at midplane. The excess in 6 = —20 to —30° is partly due to one of the 
streaks of bright areas which meet together at the South Ecliptic Pole. Thus it could be due to the so-called long term enhancement 
(Snowden et al. 1994) or scattered solar X-rays. 



yet although this problem has been known for more than 
25 years. One major reason is that there has been no en- 
ergy spectrum available in which emission line structures 
are resolved. The X-ray Imaging Spectrometer (XIS) 
(Koyama et al. 2007) on board Suzaku (Mitsuda et 
al. 2007) has a significantly improved spectral line re- 
sponse function compared to previous X-ray CCD cam- 
eras, e.g. those onboard XMM-Newton and Chandra, in 
particular below 1 keV. Although the spectral resolution 
of the instrument is not high enough to resolve the fine 
structure (triplet) of OviiKa emission, emissions from dif- 
ferent ions, e.g. Nvi, Ovii, and Oviil, can be clearly re- 
solved. Combined with the X-ray telescope (Serlemitsos 
et al. 2007), the XIS also has high sensitivity for spatially 
extended emission. Wc have observed the direction {£, b) 
= (235, 0) with Suzaku for 160 ks. The direction was se- 
lected because this point is well outside the Galactic bulge 
and north polar spur, and because the direction is an av- 
erage midplane direction without any special features, i.e. 
no bright X-sources in the XIS field of view, a typical 
neutral Hydrogen density, and a typical counting rate in 
the ROSAT all sky survey map. We found that thcOvii 
Kq, emission intensity was comparable with that of the 
MBM-12 on-cloud observation (Smith et al. 2007) and 
that a narrow bump peaked at ^ 0.9 keV was compen- 
sating the decrease of the extragalactic component. This 
strong feature, presumably due to a blend of Ne-K and 
Fc-L lines, makes the 6 = spectrum qualitatively unlike 



Table 1. Log of observation 



Observation ID 502021010 

Aim point ^ (235.00,0.00)Gaiactic 

= (113.33, -19.53)j2ooo.o 
Observation start time (UT) 20:39:20, 22 April 2007 
Observation end time (UT) 10:04:24, 25 April, 2007 
Net exposure time 189. 5ks 

^ Aim point on the focal plane was the XIS nominal po- 
sition. 

empty-field spectra at other latitudes and requires plasma 
at higher temperatures than generally seen in Galactic dif- 
fuse emission. In this paper we will show the observational 
results and discuss the origin of the excess emission. We 
will construct a model spectrum for spatially unresolved 
faint dM stars and show it can consistently explain the 
observations. 

In this paper we concentrate on the XISl data. This 
detector has a much larger effective area below 1 keV 
than the other XIS sensors because it employs backside 
illuminated CCD. Throughout this paper, we quote sin- 
gle parameter errors at the 90 % confidence level unless 
otherwise specified. 
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Fig. 2. Observed spectrum (crosses), best-fit model and its 
components (step functions) of Model (A) convolved with the 
instrument response function and residuals of the fit (bottom 
panel). The vertical error bars of data points correspond to 
the Icr statistical errors. 

2. Analysis and results 

The midplane direction, {£, b) = (235, 0), was observed 
with Suzaku in the AO-2 period. In Table 1 we show 
the log of the observation. The XIS was set to normal 
clocking mode and the data format was either 3 x 3 or 
5x5. The Spaced-raw Charge Injection (SCI) was on 
throughout the observation. We used version 2.0 pro- 
cessed Suzaku data. We first cleaned the data using the 
selection criteria, elevation from sunlit and dark earth rim 
> 20 deg, and cut off rigidity > 8 GV. We checked the 
Oxygen column density of the sunlit atmosphere in the 
line of sight of the screened data and found it to be al- 
ways below lO^^cm^^, which is the criterion for no signif- 
icant neutral O emission from Earth atmosphere (Smith 
et al. 2007). We then checked the solar wind proton flux. 
The spectrum below 1 keV could be contaminated by the 
SWCX-induced emission from the geocorona if the solar 
wind flux exceeds 4 x 10® protons s~^ cm~^ (Mitsuda et 
al. 2008). The probability of contamination is high if the 
altitude of the magnctopause is lower than ~ 10 Earth 
radii (Re) (Fujimoto et al. 2007). Here, the magnc- 
topause is defined by the lowest position along the line 
of sight whose geomagnetic field is open to interplanetary 
space. We found that the magnctopause is higher than 
10i?E throughout the observation. However we removed 
the time intervals in which the proton flux at lAU exceeds 
4 X lO'^ protons s^^ cm~^ in order to avoid any contam- 
ination by SWCX from the geocorona. After these data 
selections, the total exposure time reduced to 53 ks. 

We then constructed an X-ray image in 0.3 to 2 keV 
energy range. We detected two faint X-ray sources in the 
X-ray image and we removed circular regions centered on 
those point sources with radii of 2' and 1.5'. The radii 
were determined by the intensities of the two sources. The 
counts from the point sources outside the circular regions 
are estimated, respectively, to be less than 5 % and 3 % 
of the diffuse X-ray emission in 0.3 to 1 keV energy range. 

The non X-ray background (NXB) spectrum was con- 



structed from the dark Earth database using the standard 
method in which the cutoff rigidity distributions of the on- 
source and background data were made identical (Tawa et 
al. 2008). We found about 10 % discrepancy between the 
non X-ray background and the present observation data 
for XISl in the counting rates above 10 keV, where true 
X-ray rates are expected tobe negligible, suggesting back- 
ground uncertainty of this level. However, since in the 
energy range below 1 kcV the non X-ray background is 
only about 10 % of the diffuse X-ray background, this 
level of the background uncertainty is negligible. 

In order to perform spectral fltting, we generated an ef- 
ficiency file (arf file) for a flat field using the xissimarfgen 
software version 2007-09-22 (Ishisaki et al. 2007), as- 
suming a 20'-radius flat fleld as the input emission of the 
generator. The degradation of low energy efficiency due 
to the contamination on the XIS optical blocking filter is 
included in the arf file. The pulse height redistribution 
matrix (rmf file) was created by the script xisrmfgcn. 

In Figure 2, the energy spectrum for the accepted 
image region is shown. In the energy spectrum, one 
can clearly see a narrow bump-like structure^ peaked at 
around 0.9 keV in addition to the Ovii (0.57keV) emis- 
sion. We have carefully excluded the time intervals of the 
observation in which the data could be contaminated by 
SWCX in the upper atmosphere of the Earth (the geocoro- 
nal SWCX). However, the SWCX in interplanetary space 
(the heliospheric SWCX) cannot be excluded. The energy 
spectra of both the geocoronal and hehospheric SWCX 
emissions refiects the abundance of ions. Fujimoto et al. 
(2007) detected Fe-L and Ne-K emissions in 0.7 to 1.1 keV 
range in the geocoranal SWCX spectrum observed with 
Suzaku. Henley & Shelton (2008) suggested possibility 
of Ne-K and Mg-K emissions in SWCX based on Suzaku 
spectra. However, in both cases, emissions are detected 
as narrow lines with the energy resolution of the Suzaku 
XIS. Thus the bump structure in the present spectrum is 
qualitatively very different from those SWCX spectra. 

The narrow line features at 1.3 and 1.6 keV are instru- 
mental lines due to residuals of the NXB subtraction. We 
first fitted the spectrum in the energy range 1.1 to 5 kcV in 
order to determine the the intensities of the NXB residual 
lines. We employed a power-law function absorbed with 
a Galactic absorption for the continuum and two narrow 
Gaussian functions for the NXB lines. In the further spec- 
tral fits, we account for these instrumental features with 
two Gaussians with all parameters fixed at the values thus 
obtained. 

We then fitted the energy spectrum in the energy range 
of 0.4 to 5 keV with a model consisting of three emission 
components (Model A). In Figure 2, the best fit model 
function and its components convolved with the instru- 
ment response function are shown by step functions. The 
best fit parameter values and their statistical uncertain- 
ties are shown in the column labeled with Model (A) in 
Table 2. The first model component represents the CXB. 

^ We call this a narrow bump because this structure is broader 
than a line, but the width is only ~ 0.3 keV FWHM when it is 
approximated with a Gaussian function 
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Table 2. Best fit spectral parameters 



Component 


Model Function 


Parameter (unit) 


Model (A) 


Model (B) 


Model (C) 


CXB 


Absorption 
Power law 


Nh (10^^ cm"^) 
Photon index 


0.90 (fixed) 
1.4 (fixed) 


0.90 (fixed) 
1.4 (fixed) 


0.90 (fixed) 
1.4 (fixed) 






Normalization^ 


11 1 inn 
11. i ± U.9 


111 ^f^.,^r^A\ 
11.1 (nxeaj 


11 1 / fi,^^A \ 

11.1 (nxeaj 


LHB + 


APEC 


kT (keV) 


„ ^[.+0.051 
U.1U0_Q 


0.105 (fixed) 


0.105 (fixed) 


Heliospneric oWOX 




Normalization*^ 


1 A 1 +6.4 


lo.o ± 6.2 


11 1+2.6 


Broad line-like 

£ J. 

feature 


(V)APEC 


kT (keV) 
Normalization*^ 


766+'^-"^'-* 

U. /DD_Q 041 
„ „(.-f0.40 
' ^-0.37 


U.005_Q 072 
o on+0.45 






APEC 


abundance (Solar'^ ) 
kT (keV) 


o 1 +1.3 

0.1_i 2 


1 (fixed) 

i.OU_o.28 








Normalization*^ 




o Y+1.7 
0.1 _i2 






Absorption 


JVh (10^2 cm-2) 




0-64t°;^22 




Bremsstrahlung 


kT (keV) 
Normalization"^ 






A Q + 31.0 

'^•'^-3.7 


NXB residuals" 


Gaussian 


Centroid (keV) 


1.32 (fixed) 


1.32 (fixed) 


1.32 (fixed) 






Normalization^ 


0.24 (fixed) 


0.24 (fixed) 


0.24 (fixed) 




Gaussian 


Centroid (keV) 
Normalization" 


1.63 (fixed) 
0.23 (fixed) 


1.63 (fixed) 
0.23 (fixed) 


1.63 (fixed) 
0.23 (fixed) 


xVdof 






76.5/84 


61.6/85 


68.3/86 


^ The unit is photons 


s^^cm"^ keV'^str' 


-i@lkeV. 









^ The emission measure integrated over the line of sight ,i.e. (l/47r) J ncTiuM in the unit of 10^*cm^^ str 
" Solar abundance by Anders & Grevesse (1989). 

The emission measure integrated over the line of sight. The unit is 3.02 x lO^^^^cm^^ str^*^. 

Two intrinsically narrow Gaussians are included to represent residual instrumental emission lines. 
^ The unit is photons s~ cm~2str~*^. 



We fixed the Galactic absorption column density of this 
component to the value determined by 21 cm radio obser- 
vation (9.0 X lO^^cm-^, Dickey & Lockman 1990). The 
average AGN spectra below ~ 1 keV becomes steeper than 
above 1 keV, and the average photon index was deter- 
mined to be 1.96 below ~ 1 keV by Hasinger et al. (1993). 
We thus first tried three different models for this compo- 
nent: a power-law function with a photon index 1.4, and 
broken power-law functions with indices of either 1.54 or 
1.96 below 1.2 keV and with 1.4 above 1.2 kcV (Smith 
et al. 2007). We set the normalization of the power- 
law or the broken power-law models free. We found all 
parameters of Model (A) including the normalization of 
the CXB at 1 keV were same within the statistical errors 
for these three CXB models. Therefore we will show only 
the results with the single power-law function. The best- 
fit value of the power-law normalization is consistent with 
typical high latitude values (Revnivtsev et al. 2005). The 
intensity of the CXB component decreases rapidly below 

1.5 keV because of Galactic absorption and has negligi- 
ble contribution below 1 keV. 

The second component is a thin thermal emission 
model, APEC (http://hca-www.harvard.edu/APEC/), 
which we consider to represent the emission from hot 
gas in the local hot bubble (LHB) and from the SWCX 
process in the Heliosphere. We estimated Ovil emis- 
sion intensity by setting the O abundance of this APEC 
model to zero and substituting with a narrow Gaussian 
line at theOvii Kq, line energy, fixing the temperature of 



the emission component to the best fit value. TheOvii 
emission intensity was determined to be 2.1^o'g LU, where 
LU is photons s~^ cm"^ str~*^. This value is smaller than 
theOvii Kq emission intensity obtained from the MBM-12 
on-cloud observation by Smith et al. (2007) with Suzaku 
(3.34 ±0.26 LU). MBM-12 is located at 60 to 280 pc from 
the Earth, and theOvii emission beyond the cloud 
is completely blocked. Thus the intensity represents the 
emission from the LHB and Hcliospheric SWCX. Smith et 
al. (2007) determined the emission intensity by fitting the 
spectrum with a power-law continuum and Gaussian lines. 
However the emission intensity is dependent on how we 
describe the weak lines. Furthermore, both the data pro- 
cessing and the calibration of the instruments have been 
improved since the analysis by Smith et al. (2007) . We 
thus re-analyzed the data and determined the emission in- 
tensity using the same recipe used here. We employed the 
same models to describe the CXB and the point source 
in the field of view as Smith et al. (2007). We found 
the LHB + SWCX component of the MBM-12 on-cloud 
direction is well described with an APEC model without 
absorption with kT = 0.109!l!o'o'j'2 keV, and normalization 
= (13.4t^ j) X 10^* cm-5 str-i. The Ovil K„ emission 
was determined to be 2.93 ±0.45 LU. Both the APEC 
model parameters and the Ovii Kq emission intensities 
are consistent with the present midplane observation. 

The third component is an APEC thin thermal emission 
model at a higher temperature need to produce the bump 
peaked at around 0.9 keV. In this model the broad feature 



No. ] 

is produced by a sum of Fc-L and Nc-K lines. Wc had to 
set the the O abundance of this component free in order to 
better describe the Oviii Kq. emission. The best fit values 
of the temperature and the O abundance were respectively 
0.77 ± 0.4 keV and 3 ± 1 Solar. 

Although the value of the fit is statistically accept- 
able, we find excess emission in the 1-1.2 keV range, which 
suggests existence of even higher temperature emission. 
The large required O abundance also suggests existence 
of a lower-temperature component which emitsOviii 
emission more efficiently. Thus, the narrow bump at 
around 0.9 keV is likely multi-temperature emission. In 
the second model (Model B), we represented the narrow 
bump with two APEC components with the abundance 
fixed to 1. We found the fitting is unstable because of 
strong coupling among the two APEC components and 
the other remaining two components. We thus fixed the 
normalization of the CXB power-law function, and the 
temperature of the LHB -I- SWCX component to the best 
fit values of Model (A), respectively. In Table 2, we show 
the best fit parameters. The temperatures of the emis- 
sions were 0.66lg;Q7 and l.Slg'J keV, respectively. 

The narrow bump at around 0.9 keV may be repre- 
sented by other model functions. For example, if it is 
a sum of emissions from faint X-ray binaries, it may be 
modeled with continuum spectra. Thus, as the third 
model. Model (C), we tried to fit with a strongly absorbed 
continuum spectrum. As the continuum we adopted a 
bremsstrahlung model. As shown in Table 2, the spec- 
trum was fitted well with model parameters of kT = 
0.18^0 05 E^nd absorption column density A^h = ^ 
102icm-2. 

Finally, we divided the XIS image regions into two sub- 
regions and extracted the energy spectra separately. We 
found that the two spectra were identical to each other 
and to the total spectrum within the statistical errors. 
This suggests the emission region of the narrow bump is 
not localized within a limited image region, but spatially 
extended at least to the size of the XIS field of view (18'). 

3. Discussion 

The pointing direction of the present observation has 
no special features. Thus we can consider the spectrum 
to be representative of midplane, off-center directions. 
Therefore a similar narrow-bump emission is likely to fill 
partly the absorption of the extragalactic component in 
other midplane directions, although more observations of 
different midplane directions are necessary to confirm this 
assumption. 

First we consider spatially extended hot plasma as the 
origin. Assuming an average interstellar neutral H density 
of 1 cm~'^ and cosmic abundance, the absorption length 
of a 0.9 keV photon is estimated to be 1 kpc in midplane. 
Thus the line of depth of the observation, L, is about 
1 kpc. Then the local H density of the plasma, nn.hot, 
is estimated from the best fit model parameters of model 
(A) to be riH.hot = 1-1 x 10-^{L/1 kpc)-i/2/-i/2, where 
/ is the volume filling factor. The pressure is then p/k = 
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2.3 X 10* cm-3 K (L/1 kpc)"^^^ f"^/^ . Since / < 1, the 
pressure is likely to exceed the total midplane pressure 
derived from the vertical matter density and the vertical 
gravity, 2.2 x 10* cm'^ K (Cox 2005). We thus consider 
that a large fraction of the emission might arise from faint 
individual sources rather than diffuse gas. 

The faint sources must satisfy the following three condi- 
tions: (1) the sources have to have energy spectra which 
can be approximated either by a thin-thermal emission 
(kT 0.8 keV) or by an absorbed (TVh = 6 x lO^^cm'^) 
bremsstrahlung (kT = 0.2 keV), (2) individual sources 
must be faint and enough number of sources exist in the 
XIS field of view (18' x 18'), and (3) the sources must 
have a relatively small vertical scale height so that the 
emission is enhanced only near the Galactic plane. 

If we require at least 30 sources in the XIS field of 
view with a line-of-sight depth of 1 kpc, we need X-ray 
sources with a density higher than 0.003 pc~^. Only nor- 
mal stars and white dwarfs have high enough spatial den- 
sities. Among those objects, the contribution of main se- 
quence M stars (dM stars) dominates the emission in the 
0.3 to 1 keV energy range because they have relatively 
high X-ray luminosities and high number density (Kuntz 
& Snowden 2001). Young (^ 1 Gy) stars have a large 
X-ray luminosity and dominate the total X-ray emission. 
They have a small vertical scale heights (a few lO's pc, 
e.g. Bienayme, Robin & Creze 1987). 

Giampapa et al. (1996) determined the energy spec- 
tra of six nearby (< 10 pc) dM stars with ROSAT PSPC 
observations. The spectra were found to consist of two 
thermal components. The average temperatures were 
kT = 0.138 and 0.78 keV, respectively, when they were fit- 
ted with solar-abundance plasma models. The two compo- 
nents have comparable emission measures. The high reso- 
lution spectrometers onboard Chandra and XMM-Newton 
observatories resolved the emission spectra into number 
of emission lines (Maggio et al. 2004; Giidel et al. 
2004). Maggio et al. (2004) determined the emission 
measure distribution (EMD) as a function of temperature 
for the dM3e star, AD Leo. The EMD show two peaks at 
logT = 6.2 and logT ^ 6.9, thus at kT = 0.14 keV and 
0.68 keV, which is consistent with the two temperature 
model from the ROSAT observations. The temperature of 
the higher-temperature component is consistent with that 
producing the narrow bump in the spectrum of present ob- 
servation. The lower temperature component of dM stars 
emits mostly Ovii Kq in the XIS band. Because opacity 
at Ovii is about twice as high as that for 0.9 keV photons, 
the line-of-sight depth contributing to the present obser- 
vation is a factor of two shorter. Thus if we assume the 
narrow-bump emission arises from the 0.78 keV compo- 
nent of dM stars, we find the Ovii Kq, emission intensity 
from the 0.14 keV component of dM star is only about 
30% of that of the LHB + HS-SWCX component, which 
is within the statistical error of the intensity. Thus the 
existence of the M-star's lower-temperature component in 
the present observed spectrum is statistically acceptable. 

Giampapa et al. (1996) also suggested that the inten- 
sity of the higher temperature component shows positive 



Unresolved Soft X-ray Emission from the Galactic Disk 



6 



Masui et al. 



0.1 



.4:-- ! : 
" L.l ■ -.1""' 


1 1 1 1 

■ : dM stars •] y 
"(i 30% of model) 1 : : .- 1; 

^ ! '■■ 

r ! 1 ! : 
^1 .... 1 
■Hi ■ ' I ■ 






- 1 




m 





0.5 1 2 

Energy (keV) 

Fig. 3. Observed spectrum (crosses), best-fit model and its 
components (step functions) convolved with the instrument 
response function and residuals of the fit (bottom panel). 
Here the high-temperature component of Model (A) was re- 
placed with the model spectrum of faint dM stars. The verti- 
cal error bars of data points correspond to the Itr statistical 
errors. 
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Fig. 4. Residuals of R45 band surface brightness from the 
unabsorbed constant emission + CXB model in Figure 1 (step 
functions), and expected X-ray fluxes from faint dM stars 
(solid curve) as functions of the Galactic latitude, along 
£ = 235°. The surface brightness of the dM star component 
inferred from the spectral flt to the Suzaku spectrum is shown 
as the fiUed-circle data point at 6 = 0. 



correlation with X-ray intensity during flare events. The 
EMD of AD leo obtained with the Chandra grating ob- 
servation is extended at temperatures above the peak at 
0.68 keV up to - 2 keV. Van den Bcssclaar et al. (2003) 
showed that the EMD of AD Leo is significantly enhanced 
in the temperature range of kT = 1-2 keV during the de- 
cay of a flare. Wargelin et al. (2008) detected a flare 
from the nearby M star, Ross 154 with Chandra ACIS 
and found that the temperature of the higher temperature 
component increased during the flare when it was fitted 
with a two component model. During the decay phase of 
the flare, the temperature was 1.9 keV and the emission 
measure increased by a factor of about three compared to 
that during quiescence. The higher temperature (kT = 
1.5 keV) component in Model (B) has an emission mea- 
sure comparable to that of 0.7 keV component. This could 
be the contribution of high temperature emission in EMD 
of dM stars in quiescence and flare. The emission measure 
in 1 to 2 keV range is about 0.3 times that in 0.7 to 0.9 
keV range for AD Leo. Then if about 20 % of dM stars 
in the field of view are in flaring state, the 1.5 keV com- 
ponent can be explained. Since the dM stars are highly 
variable (e.g. Collura et al. 1988, Schmitt 1994) , we 
consider this is plausible. 

The contributions of stars to the soft X-ray back- 
ground have been estimated by several authors (Rosncr 
et al. 1981; Schmitt & Snowden 1990; Guillout et al. 
1996; Kuntz & Snowden 2001). However, they estimated 
contributions only for high latitudes (6 > 20°). The X-ray 
luminosities of stars are distributed in a relatively narrow 
range and the distribution function can be approximated 
by a log-normal function, when they are sorted by spectral 
type and age (Schmitt & Snowden 1990; Schmitt et al. 
1995) . Assuming a double-component thin thermal emis- 
sion of kT = 0.766 keV (the best fit value of Model A) and 
0.138 keV with the same emission measures for the two 



components, we can convert the the average luminosities 
of the log-normal distributions to the average emission 
measures. By multiplying them with the midplane star 
densities and line of sight length contributing to the dif- 
fuse emission (^ 1 kpc), and adding them together, we can 
roughly estimate the total emission measure. Using the 
values compiled by Kuntz & Snowden (2001) for luminosi- 
ties and star densities, we obtained 4.1 x lO^"' cm^^ str^^ 
(Table 3). This is consistent with the observed value of the 
emission measure, (3.8 ±0.4) x 10^"* cm~^ str~^ (Table 2). 
Among different age subgroups, early and late type dM 
stars with ages of 0-1.5 Gy contribute a large fraction 
(^70 %) of the emission measure. The average number 
of sources of this age group is 60 in the XIS field view. 

The energy flux of the faintest class of the dM stars 
is estimated to be 6 x 10~^^ erg s~^ cm~^, assuming 
the average luminosity of the age class, a 1 kpc dis- 
tance, and the average Galactic absorption. The sensitiv- 
ity of the Chandra Multiwavelength Plane (ChaMPlane) 
Survey reaches close to the limit. However, since the 
ChaMPlane fields contain directions with Galactic lati- 
tude up tp \b\ = 10°, the logA^- logs' relation in the 0.4-2 
keV energy band is significantly contaminated by extra- 
galactic sources (Hong et al. 2005). Moreover, as we 
will show later the surface density of the dM star drops 
rapidly with |&|. We thus estimate the number of sources 
in the XIS field using the source counts in 0.4 - 2 keV 
range of the two low-latitude ChaMPlane fields: ObsID 
2810 {b = 0.18°) and ObsID 2218 {b = 0.14°) which have 
modest exposure times: 48.8 ks and 30.2 ks, respectively. 
Simply scaling the source counts of the level 3 analysis 
of Hong et al. (2005) by the ratio of the field of views, 
wc obtained 70 and 24 as the expected source counts for 
the XIS field of view. The nature of those sources are yet 
not known. However these source numbers are consistent 
with the above estimation of dM stars, i.e. 60. 
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Table 3. Rough estimation of midplane omission measures of dM stars 



Type 


Age 


< logLx 


cr(logLx 


a 


<EM>'^ 


Density"^ 


Total EM^ 




(Gy) 


(log erg s"^) 


(log erg s~ 


') 


(lO^^cm-^) 


(10"3pc-3) 


(lO^'^cm-^str-i) 


M (early) 


0-0.15 


29.19 


0.32 




48.3 


4.28 


1.72 




0.15-1 


27.89 


0.72 




7.30 


8.96 


0.55 




1-10 


26.86 


0.77 




0.830 


3.29 


0.23 


M(late) 


0-0.15 


29.12 


0.49 




59.3 


2.45 


1.21 




0.15-1 


28.2 


0.49 




7.12 


5.12 


0.30 




1-10 


27.22 


0.49 




0.746 


18.8 


0.12 



Sum 4.13 

^ Log normal probability distribution function of M star luminosity (Kuntz & Snowden 2001). 
^ Emission measure for the average luminosity of the log-normal distribution. Double com- 
ponents of the same emission measure with temperatures of kT = 0.138 and 0.766keV are 
assumed. 

Average midplane star density (Kuntz & Snowden 2001). 
Total emission measure along a line of sight of 1 kpc length. 



We then constructed model energy spectra, taking into 
account the Galactic absorption and vertical star distri- 
bution, in order to compare the emission spectrum with 
observation, and to estimate the dependence of the surface 
brightness on Galactic latitude. For the neutral hydrogen 
column density, we used Equations (1) (molecular) and 
(2) (atomic) of Ferriere (2001). We scaled the midplane 
densities at the various galactocentric radii by the surface 
density in Figure 1 of the same paper. For molecular 
gas we used the curve from Clemens et al. (1988) in the 
figure. The atomic gas surface density depends on galac- 
tocentric radius only outside 14 kpc; thus it does not af- 
fect the model spectrum constructed below. For the star 
density, we follow Kuntz & Snowden (2001), and used 
Equations (Al) (0 - 0.15 Gyr old stars) and (A2) (0.15 
- 1 Gyr old stars) of Bicnaymc, Robin & Crezc (1987). 
Then the emission spectrum is given by 

1 f°° 

f{E) = — / dsy{A{E,TH) + A{E,n))EM,p,{s) 
Jd„ V 

xexp -crabs(-E) / ds'iiuis') 

where the integrations with s and s' are done along the 
line of sight, while the summation by i is taken for dif- 
ferent classes, i.e. spectral types and age groups, of dM 
stars. EMi and pi are, respectively, the emission measure 
of the average luminosity star and spatial density of class 
i. The emission spectrum from a plasma with a tempera- 
ture of T and unit emission measure is given by A(_E,T), 
and (Tabs and nn are respectively the X-ray absorption co- 
efficient and the neutral hydrogen density. Finally, do is 
the distance within which dM stars are resolved as indi- 
vidual stars. We varied it from 5 pc to 20 pc and found 
there was no significant difference in the resultant model 
spectra. 

We performed spectral fits using the model function 
f{E) instead of a thin thermal emission model for the nar- 
row bump in Model (A). In the fit, we introduced a scaling 
parameter for f{E) and set it free. We also set the normal- 



ization factors of the LHB+SWCX and the CXB compo- 
nents free. The result is shown in Figure 3. A minimum 

value of 87.87 with 87 degrees of freedom was obtained 
for a scaling parameter value of 1.27j;Q |^2- The normal- 
ization of both the CXB component (10.1^Q g'°) and the 
LHB+SWCX component (13.2 ± 3.0) were smaller than 
that of Model (A), however their differences are within 
the statistical errors. Therefore, both the surface bright- 
ness and spectrum of the diffuse emission at {£,b) — (235, 
0) can be consistently explained by introducing the unre- 
solved emission from dM stars. 

In Figure 4, we show the dM star model fluxes in the 
ROSAT R45 band as functions of 6 for ^ = 235°. We also 
show the residuals of the observed ROSAT R45 band flux 
from the unabsorbed constant emission + CXB model in 
Figure 1. As shown in this figure, the dM star emissision 
can fill in the notch in |6| < 2°, as well as providing the 
high-temperature feature seen in the Suzaku observation. 
However, there exists more excess flux at b^2 — 10° than 
the dM star model predicts. Thus yet another emission 
component which should have a double-peaked latitude 
structure must fill the notch. More Suzaku observations 
in |6| = 2° — 10° are necessary to elucidate this unknown 
emission component. 
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